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ABSTRACT 

m ■ „ 

Context. This is the first step in a project to combine studies of eclipsing binaries and oscillating stars to probe the interior of Blue Stragglers 
f-^. (BS). This may imply a way to discriminate observationally between different birth mechanisms of BS stars. 

f^**) . Aims. We study the open cluster NGC2506 which contains oscillating BS stars and detached eclipsing binaries for which accurate parameters 
can be derived. This will tightly constrain the cluster isochrone and provide an absolute mass, radius and luminosity-scale for the cluster stars 
along with the cluster age, metallicity and distance. The present work focuses on obtaining the light curves of the binaries and determine their 
I • orbital periods, on obtaining power spectra of the oscillating BS stars to select targets for follow-up studies, and on searching for y Doradus 
O ' type variables which are also expected to be present in the cluster. 

+3 , Methods. With a two-colour, dual-site photometric campaign we obtained 3120 CCD-images of NGC2506 spread over four months. We 
^2 analysed the BI time-series of the oscillating stars and used simulations to derive statistical uncertainties of the resulting frequencies, amplitudes 
and phases. A preliminary mode-identification was performed using frequency ratios for the oscillating BS stars, and amplitude ratios and phase 
^ . differences for a population of newly detected y Doradus stars. 
J^j Results. We quadrupled the number of known variables in NGC2506 by discovering 3 new oscillating BS stars, 15 y Doradus stars and four 
rS , new eclipsing binaries. The orbital periods of 2 known, detached eclipsing binaries were derived. We discovered a BS star with both p-mode 

■ and g-mode variability and we confronted our y Doradus observations with state-of-the-art seismic models, but found significant discrepancy 
between theory and observations. 

Conclusions. NGC2506 is an excellent target for asteroseismic tests of stellar models, as strong external constraints can be imposed on the 
models of a population of more than 20 oscillating stars of different types. 

Key words, stars: blue Stragglers, stars: variables: delta Set, stars: binaries: eclipsing - open clusters and associations: individual: NGC2506 - 
techniques: photometric 



1. Introduction 

Blue straggler (BS) stars are still a much debated phenomenon. 
In a cluster HR diagram, these stars are situated above and blue- 
ward of the main sequence (MS) turn-off point, where more 
massive stars were located before they evolved away to red gi- 
ants. Their appearance in a region of the HR diagram where 
they should not be contradicts our standard single-star evolu- 
tion theory. 
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Although the origin of BS stars is still not fully understood, 
the generally accepted working hypothesis is that these stars 
originate from stellar collisions, binary coalescence, or from 
mass transfer in a binary system. These ideas go back to the 
pioneering work of Hills & Day (1976), and McCrea (1964). 
Because of these stellar interactions, a BS star is formed with 
a mass greater than the turn-off mass of the cluster. Such a star 
is therefore able to prolong its life on the MS, or even reappear 
as a born-again star, which explains why BS stars have not yet 
evolved far away from the MS. 

However, before this hypothesis can be promoted to a the- 
ory, more observational and theoretical evidence needs to sup- 
port and quantify it, and this is what current BS research is fo- 
cusing on. We mention the recent result of Ferraro et al. (2004) 
who found a bimodal radial BS population distribution in the 
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globular cluster 47 Tuc, which they interpreted as evidence for 
different but simultaneously working BS birth mechanisms (see 
also Mapelli et al. 2006 for a follow-up article). For the same 
cluster Ferraro et al. (2006) pointed out that at least some of the 
BS stars seem to show the chemical signature of binary mass 
transfer, adding another piece of evidence for the mass trans- 
fer scenario. Progress on understanding the collisional scenario 
has mainly been made on the theoretical side. Sills et al. (2001, 
2005) used smoothed particle hydrodynamics simulations to 
show that off-axis collision products have rotational velocities 
well above their break-up velocity so that they can only survive 
if they can somehow get rid of the excess angular momentum. 
The authors mention magnetic breaking caused by a circum- 
stellar disk as a possible cause of a BS spin down. Such discs 
may (De Marco et al. 2004) or may not (Porter & Townsend 
2005) already have been discovered around BS stars. 

Clearly, the theories for the different BS birth mechanisms 
are getting more established through the increasing amount of 
evidence. It's still difficult, though, to recognize how a partic- 
ular BS star was born. Especially a diagnostic to discriminate 
observationally between binary coalescence and stellar colli- 
sion is not available at the moment. This makes it difficult to 
study, for example, to what extent and for how long the inte- 
rior of a BS star formed by "unrolling" a donor star is different 
from one formed by a direct collision of two stars. 

Might asteroseismology help in this respect? In many clus- 
ters the BS stars are situated in the classical instability strip 
where one expects 6 Scuti or SX Phe oscillations. Many oscil- 
lating BS stars have indeed already been detected, see for ex- 
ample Gilliland & Brown (1992), Gilliland et al. (1998). These 
oscillations could provide a way to model the BS star interior. 
The fact that the BS stars are part of a cluster makes them even 
more suitable asteroseismic targets, as chemical composition 
and age may be assumed constant for each BS member. These 
constraints can be made quantitatively more stringent if the 
cluster also contains detached eclipsing binaries for which ac- 
curate masses and radii can be determined, so that the isochrone 
can be precisely calibrated (see Grundahl et al. 2006). Finally, 
a combined set of observed frequencies in a number of oscil- 
lating stars within the same cluster, can be used to extract more 
information than for single stars (Frandsen & Kjeldsen 1993). 

This paper reports on the first stage of this Blue Straggler 
project: the challenging task of deriving useful oscillation spec- 
tra of BS stars in the open cluster NGC2506. The reason for 
choosing NGC2506 is fourfold. First, the cluster contains BS 
stars of which three of them were shown to oscillate by Kim 
et al. (2001), by means of photometric time series (V1-V3). 
Secondly, the cluster contains detached eclipsing binaries. One 
was found by Kim et al. (2001; V4), and another by ourselves 
(V5). Thirdly, it is an open cluster close to the celestial equator 
so that it is visible from both hemispheres, making a multi-site 
campaign of variable stars possible. Finally, most of the cluster 
is confined within a radius of 6 arcmin, making it fit nicely in 
a typical CCD image; yet the cluster is not so crowded that the 
photometric precision is affected. 

Besides the BS project, we present another project which 
started as a side result but which has meanwhile grown into 
an equally important project. The turn-off stars of NGC 2506 



span a large part of the y Doradus instability strip making the 
cluster perfect for investigating the y Doradus phenomenon, 
should such stars be present. Kim et al. (2001) searched for 
y Doradus stars without detecting any. The data presented here 
are both more extensive and have a higher precision per data 
point, making the search more sensitive and, as it turns out, 
more successful because we do detect a surprisingly high num- 
ber of y Doradus stars in the cluster. These y Doradus stars 
are very promising asteroseismic objects as their high-order 
g-modes probe deeply in the stellar interior. In addition, they 
are exciting laboratories to study stellar physics such as semi- 
convection, convective overshooting, and (both slow and fast) 
stellar rotation. 

Only recently, however, could some quantitative compar- 
ison between theory and observations be carried out for the 
y Doradus stars. Dupret et al. (2004, 2005a) explained the 
driving mechanism of the y Doradus g-modes by using the 
time-dependent convection (TDC) treatment of Grigahcene et 
al. (2005). With the same TDC models, Dupret et al. (2005b) 
explained for the first time and with success the photometric 
amplitude ratios and phase differences of five field y Doradus 
stars. In the same way as for the BS stars, our cluster y Doradus 
stars provide a more stringent test of the theory. This paper re- 
ports the discovery of many new y Doradus stars, the observed 
frequencies, and a first attempt to model the photometric am- 
plitude ratios and phase differences with the TDC theory, with 
remarkable results. 



2. The target cluster: NGC2506 

NGC2506 (CK2000, £2000 = 08 h 00 m 01 s , -10°46'l2") is an old 
open cluster that already received quite some attention in the 
literature. The morphology study of Chen et al. (2004) indi- 
cates that it is a mildly elongated cluster counting about 1091 
stars, and the recently published catalogue of Kharchenko et 
al. (2005) puts NGC2506 at a distance of 3460 pc. There 
seems to be a consensus in the literature about the age of the 
cluster, which would be t = 2.14 + 0.35 Gyr (Salaris et al. 
2004) although Xin & Deng (2005) quotes an age of 3.4 Gyr. 
Such an old age is inconsistent with our Stromgren photom- 
etry which also points to an age close to 2 Gyr. We refer to 
Carretta et al. (2004) for a literature overview on the metallic- 
ity of NGC2506. Although the different values vary quite a lot, 
all chemical analyses show that the cluster is metal deficient. 
We give a somewhat higher weight to the value obtained by 
Carretta et al., [Fe/H] = -0.20 + 0.01, as it is the only metal- 
licity obtained through an iron abundance analysis using high- 
resolution spectroscopy. Nevertheless, we treat their quoted er- 
ror bar as a lower limit, as their average metallicity value is 
based on a sample of only 2 stars. Also, their value differs from 
other values reported in the literature obtained from photome- 
try or low-resolution spectroscopy. Gratton (2000), for exam- 
ple, lists a value of [Fe/H] = -0.42 + 0.09. This leads us to 
believe that the uncertainty on the metallicity is actually larger 
than 0. 1 dex. As a final cluster parameter, we adopt the redden- 
ing given by Carretta et al. (2004): E(B-V) = 0.073 + 0.009. 



T. Arentoft et al.: Oscillating Blue Stragglers, y Dor stars and EBs in NGC2506 



3 



3. The data and the data reduction 

Our photometric time-series CCD observations were obtained 
using DFOSC at the Danish 1.54-m telescope at ESO, La Silla, 
and MEROPE at the Flemish 1 .2-m Mercator telescope at La 
Palma. At ESO we observed during 23 nights in January 2005, 
at La Palma during 20 nights in the period January- April 2005. 
Overlapping data were obtained on some of the January nights. 

In order to optimise the time series to the highest possible 
precision, care was taken to maintain the same pointing during 
the observations, i.e., to keep the stars at fixed pixel positions 
on the CCD. This was done by starting each night with the same 
pointing as the previous nights (within a few pixels) and keep- 
ing this pointing using autoguiding. The exposure times were 
chosen so that the brightest of the known variables, VI, were 
below the saturation limit. At ESO, this led to exposure times 
of approximately 100 sec in B and 30 sec in /, hence sequences 
of (2B, 41, 2B, ...) were obtained. At the Mercator, longer expo- 
sure times were possible due to the smaller telescope aperture, 
but 120 sec were used for both filters in order to maintain suf- 
ficient time resolution for the short-period oscillating BS stars. 

In total, 2100 BI images were obtained at ESO, and 1020 
at La Palma, corresponding to 101.8 (B) and 97.8 (/) hours of 
observing time at ESO, and 60.9 (B) and 60.4 (I) at La Palma. 
In Fig.[T]all data are shown for the variable V3, illustrating the 
time distribution of the data set. Most of the data were obtained 
in January 2005 (136/132 hours for B/I) with additional data 
obtained in March and April at the Mercator. In the gap be- 
tween the first and second batch of observations (Fig.Q] upper 
panels) a few data points were obtained per night at La Silla, to 
monitor the EBs. These data are not included in the analysis of 
the oscillating stars presented here. The data from March and 
April were obtained at higher airmass and had lower exposure 
levels (#counts/star), resulting in higher noise per data point. 

The CCD images were calibrated using standard proce- 
dures. The mean BIAS level was subtracted from each frame 
using overscan areas, and a median BIAS image (with zero 
mean) was then subtracted. Evening and morning sky-flats had 
been obtained whenever possible, and from those, a single, me- 
dian sky-flat were constructed per filter and per observing run. 
The run at the D1.5-m spans nearly a month, but before apply- 
ing the flat-field calibration to the science images, we checked 
that the flat-fields were indeed stable over the length of the ob- 
serving run. 

The photometric reductions were done separately for the 
Mercator and ESO data, using the software package MOMF 
(Kjeldsen & Frandsen 1992). MOMF applies a very robust al- 
gorithm combining PSF and aperture photometry, and is ex- 
cellent for time-series photometry in semi-crowded fields. The 
data from each site were merged on a star-by-star basis, by 
shifting the mean levels of the individual light curves to zero. 
An example of this merging is shown in Fig. [2] for V3 dur- 
ing a single night. Not all stars were present in both data sets; 
due to differences in field-of-view (MEROPEs FOV is 6.5x6.5, 
DFOSCs 15'xl5') more stars were observed at ESO. The ob- 
served field is shown in Fig. [20] (a DFOSC image), which also 
serves as a finding chart for the variable stars discussed below. 
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Fig. 1. The total data set for the variable V3, showing the time 
distribution of the combined data from the two sites (ESO and 
La Palma). 

3.1. Analysis of differential light curves 

The output from MOMF is differential time series photometry 
for all stars in the images, using the mean of all stars as the 
reference level to take out atmospheric effects. The light curves 
were kept in the instrumental system, but calibrated Stromgren 
photometry has previously been obtained for the cluster by F. 
Grundahl; these data are presently unpublished, but calibrated 
as described in Grundahl et al. (2002). We will in the present 
paper rely on these data for deriving stellar parameters for the 
variable stars. 

We obtained light curves for 863 stars in the MEROPE field 
and 3354 stars in the DFOSC field, of which one must expect 
that a significant fraction, especially of those far from the clus- 
ter center, are not physical members of the cluster. The result- 
ing photometric precision in the light curves was as aimed for, 
with a precision per data point of 1-2 mmag for the bright stars 
(on good nights), gradually decreasing for stars of fainter mag- 
nitudes as the photon noise increases. 

The aim of our data analysis was to detect all variable stars 
in the observed fields, and to use all data for a detailed analysis 
of the variability. The first step in the analysis was to clean all 
light curves for obvious outliers and to cross-identify the com- 
mon stars in the MEROPE and DFOSC fields. We then cal- 
culated amplitude spectra for all stars, based on the MEROPE 
and DFOSC data separately and inspected the light curves and 
amplitude spectra of each star in the two data sets manually, to 
select possible variable stars for further analysis. For stars to be 
selected as candidates, we required that variability was present 
in both the MEROPE and DFOSC data (for common stars) and 
in both the B- and / filter data. 

At this stage we had a list of 33 variable star candidates, 
of which 28 ended up being retained as variables, while 5 stars 
were rejected at a later stage in the analysis. For the candidates 
observed at both telescopes we now merged the MEROPE and 
DFOSC data by shifting the light curves to zero mean. For 
some of the blue stragglers (such as V3; see Fig.|2|, we sub- 
tracted some incoherent, low-frequency variations in order to 
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Fig. 2. Overlapping data on V3. Triangles are data from La 
Palma, squares are from La Silla. A slow trend with a frequency 
of 1.0 d -1 was present in the La Silla B-data only, and was sub- 
tracted before merging the data (see Sections 3.1 and 4.1). 



match up the data from the two sites. As these stars are bluer 
than the large majority of cluster stars, we ascribe these varia- 
tions to second-order extinction effects that remain uncorrected 
because the differential photometry is dominated by the redder 
cluster turn-off stars. 

Because the data quality during a time series is variable, 
we assigned statistical weights to each data point (see e.g., 
Arentoft et al. 1998, Handler 2003). This was done by highpass 
filtering the time series, using the Fourier transform, so that 
the filtered light curve only contained high-frequency (white) 
noise, which is useful for finding deviating points. We tried 
different weighting schemes based on the deviation of the in- 
dividual data points from the light curve mean compared to 
the overall standard deviation (cr) of the light curve. In the end 
we chose to give points deviating less than 2cr in the highpass 
filtered series weight 1.0, while points deviating more were 
weighted as l/o -2 , as this scheme resulted in the lowest noise 
in the weighted amplitude spectra. 

We then scrutinized the evidence for variability, using 
weighted amplitude spectra as well as the light curves, for 
each of our 33 candidates. In this way we detected 28 vari- 
ables, of which only 5 were previously known. The variables, 
which will be discussed in detail below, are 6 oscillating blue 
straggler stars (6 Scuti oscillations), 15 y Doradus candidates, 
6 binaries and one star of presently unknown type. The clas- 
sification is based on the time scale of the variability (~days 
for the y Doradus candidates, ~hours for the oscillating blue 
stragglers), BI amplitude ratios - to separate y Doradus stars 
from, e.g., ellipsoidal variables - as well as the position of the 
stars in the colour-magnitude diagram. Fig. [3] plots this dia- 
gram based on the Stromgren data mentioned above, with the 
detected variables and the observational 6 Scuti (Breger 2000) 
and y Doradus (Handler & Shobbrook 2002) instability strips 
indicated. In fact, we will in the following argue, that the ma- 
jority of the y Doradus candidates, which are all positioned in- 



side the instability strip, can be considered bona fide y Doradus 
stars (Handler 1999, Handler & Shobbrook 2002). 

The last step in the data analysis was a detailed frequency 
analysis of the oscillating stars, and a period search for the 
eclipsing binaries. For the oscillating stars, the frequency anal- 
ysis was carried out using Period98 (Sperl 1998), which uses 
simultaneous least-squares fits of all detected frequencies in 
a light curve. We only include frequencies that are present in 
both the B and the / data, with a Signal-to-Noise ratio in the 
amplitude spectra of 4 in at least one of the filters, although 
most detected frequencies meet this requirement for both fil- 
ters. Because our B and / data have different sampling, the 
spectral window functions are slightly different for the two 
cases. This may lead to small differences in the detected fre- 
quency values between the two filters. For each detected fre- 
quency, we therefore manually inspected the amplitude spectra, 
to find a frequency value that described the variation in both fil- 
ters well. For each star, we then kept the detected frequencies 
fixed, and fitted all frequencies simultanously to the B and I- 
band data separately, to determine amplitudes and phases for 
both filters. 

All detected frequencies for the oscillating blue stragglers 
and the y Doradus stars are given in Tables Q] and [2] and 
Stromgren indices for all variables, as well as stellar parameters 
for most of these stars, have been collected in Table[3] We have 
used the calibration of Alonso et al. (1996) for deriving effec- 
tive temperatures, the bolometric corrections follow Lejeune et 
al. (1998), and luminosities have been derived using the known 
cluster distance and Mb i,o = 4.746. 

3.2. Mode parameter uncertainties 

As we plan to use the amplitude ratios and phase differences for 
theoretical modeling of the y Doradus stars, we should have 
an idea how large the uncertainty is on the observed values. 
One possibility would be to use the well-known expressions 
for the standard errors of the amplitudes and phases given by 
e.g. Breger et al. (1999). One should be aware, however, that 
these expressions were derived for equidistant time series, and 
under the assumption that the true oscillation frequencies are 
known. Our time series are notoriously non-equidistant. In ad- 
dition, as many of our targets are multi-periodic, we can expect 
that the derived frequencies for the lowest-amplitude modes 
may deviate somewhat from the true eigenfrequencies. Hence, 
the standard formulas for calculating errors are not directly ap- 
plicable to our data set. 

We therefore used extensive simulations to obtain uncer- 
tainty estimates. For each star we simulated B and / lightcurves 
with exactly the same time sampling as for the observed time 
series. Each time series consisted of a sum of sines with exactly 
the same frequencies, amplitudes and phases as derived from 
the observed time series. Next, white noise was added with the 
same rms value as in the residuals of the observed time series. 
These time series were then subjected to an automated analysis 
procedure that mimicked accurately how we analysed the ob- 
served time series. Oscillation frequencies were searched for 
in the B time series with cyclic prewhitening, and these fre- 
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Fig. 3. Colour-magnitude diagram for NGC2506 with the de- 
tected variables marked. The blue dashed lines are the observa- 
tional 5 Scuti instability strip (Breger 2000) and the red solid 
lines the y Doradus instability strip from Handler & Shobbrook 
(2002). A distance modulus of 12.65 and E(b -y) = 0.054 was 
used for placing the instability strips in the CMD. Blue squares 
marks the oscillating BS stars, red triangles the y Doradus can- 
didates, and black stars the eclipsing binaries. The insert in the 
lower left corner shows a zoomed view of the y Dor instability 
strip. (Colours are only visible in the electronic version). 

quency values were then used to perform a least-squares fit of 
the B and / time series to obtain the amplitudes and phases. The 
frequency analysis could be automated because the algorithm 
searched in predefined intervals [v; - j, V; + j] where v, is the i- 
th observed frequency (used as input to construct the synthetic 
time series) and T is the total time span of the time series. From 
the result of this fit, amplitude ratios and phase difference could 
be computed. For each star, the procedure was then repeated 
4000 times. The result is for each oscillation mode of each star, 
a set of 4000 amplitude ratios and a set of 4000 phase differ- 
ences, one for each noise realisation in the time domain. Note 
that because of the setup with the predefined intervals, the au- 
tomated procedure cannot be fooled by one-day aliases. It can 
still be misled, however, by alias frequencies smaller than 2/T. 

There are two questions we would like to answer with these 
simulations. First, can we retrieve the frequencies, amplitudes 
and phase that were used as input to construct the synthetic time 
series? This is a necessary condition to have confidence in the 
observational results, which are listed in Tables Q] and [2] The 
short answer is "yes". The sets of simulations can be used to 
construct histograms of the fit parameters, which in turn can be 
used to compute average values and standard deviations. The 
a posteriori derived mean frequencies are always very close to 
the input frequencies, which is reassuring. Sometimes, how- 
ever, the histogram is bimodal or even trimodal with well sep- 
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Fig. 4. The amplitude ratio I/B versus the phase difference I—B 
for all oscillation frequencies detected in the y Doradus stars. 
The coordinates of the bullets are the observed values, the error 
bars were computed from the simulated time series. Red trian- 
gles are used for those highest-amplitude modes for which no 
aliasing is observed in the simulated time series. 

arated peaks. For these frequencies the problem of aliasing is 
very severe, and this should be taken into account when deter- 
mining the uncertainties. 

Second question: how large are the uncertainties of the ob- 
served amplitude ratios and phase differences? We use the vari- 
ances derived from the simulated datasets as an uncertainty 
measure for the observed quantities. Always the total variance, 
i.e. intra-variance + inter-variance, was computed taking into 
account possible confusion with window aliases. The results 
are shown in Fig. [4] for the y Doradus stars. For every target we 
checked whether the highest-amplitude mode showed aliasing 
in the simulated data (i.e. a multi-modal histogram) and desig- 
nated the ones that don't, with red triangles in the figure. We 
regard the observational amplitude ratios and phase differences 
of these latter modes as most reliable. As we expected, all our 
uncertainty estimates are comparable or somewhat larger than 
the estimates we would have obtained with the formalism of 
Breger et al. (1999). We postpone futher discussion of Fig.|4]to 
Section 5. 

As a final comment we should mention that for 3 modes (vi 
of V12, V2 of V15, and v\ of V17) the phase differences showed 
a weak trend with respect to the frequencies in the dataset of 
4000 realisations. This signals that we should treat the obser- 
vational values for these 3 modes with caution. 

4. Oscillating blue straggler stars 

NGC2506 has a known population of 12 BS stars, which are 
all probable cluster members based on proper motions (Xin 
& Deng 2005, Chiu & van Altena 1981). Kim et al. (2001) 
found 3 oscillating BS stars (V1-V3), of which one (V3), 
was not included in the proper motion study of Chiu & van 
Altena. We discover an additional three oscillating BS stars 
(V6-V8) of which V7 has a membership probability of 42% 
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from the proper motion study, while the proper motions of 
the other two were not measured. Although we do not have 
spectroscopic data to pinpoint their membership status, all six 
stars have Stromgren indices in agreement with cluster mem- 
bership - their m\ indices, for instance, are similar to those 
of the cluster turn-off stars. The positions of the six stars in 
the CMD are marked with blue squares in Fig. [3] and the ba- 
sic parameters are given in Table [3] Furthermore, the six stars 
have short oscillation periods (Table Q] frequencies around 10- 
20 d ), are mostly multiperiodic and, except V2, have low 
amplitudes (mmag). These are the signatures of 5 Scuti type 
oscillations. Because the 6 Scuti instability strip overlaps the 
region of the blue stragglers (Fig. [3]) the six stars are most likely 
cluster members, as an instrinsically fainter foreground star, or 
a brighter background star, would not be expected to display 
6 Scuti type oscillations. 

Light curves from a single night for all 6 stars are shown in 
Fig. [5] V6 and V8 are outside the MEROPE FOV so we have 
only single site data for these stars. This results in a different 
spectral window function, however for both the single- and the 
dual-site data, prominent sidelobes are present in the ampli- 
tude spectra. We show such spectra of V1-V3 in Fig.[6]and of 
V6-V8 in Fig. [7] The window function for the stars observed at 
both sites is best illustrated by V2, which has a single dominant 
frequency, while that of the single-site data is illustrated by V6. 
For this star, the signature of the dominant frequency is, how- 
ever, distorted by a low-amplitude mode at higher frequency. 
Below, we discuss the individual stars in more detail. 

4.1. The known variables V1-V3 

Kim et al. (2001) detected short-period oscillations in V1-V3, 
with dominant frequencies of 14.742, 10.854 and 12.264 d _I , 
respectively, in good agreement with our results. They clas- 
sified all three variables as 6 Scuti stars, but only VI as a 
blue straggler. However, based on our new, more precise CMD 
(Fig- Ell, the six short-period variables are bluer than the cluster 
turn-off and should therefore all be classified as blue stragglers. 
Using a period-luminosity relation for 5 Scuti stars, Kim et al. 
arrived at a tentative identification for the dominant modes of 
VI and V2 as 3rd radial overtone and the fundamental radial 
mode, respectively. 

VI is the brightest and bluest star of the variables. The 
light curves of this star were corrected for low-frequency (0- 
2 d~') variations at the level of 2-3 mmag before the data were 
merged. These variations caused minor mis-matches between 
the MEROPE and DFOSC data, and must therefore be ascribed 
to instrumental or atmospheric effects (see Sect. 3.1). In the fi- 
nal light curves, we detect 7 frequencies in the range 13—18 d _1 , 
with amplitudes between 2 and 5 mmag (TableHJ. V2 is a high- 
amplitude variable (0. 12 mag in B). The data from the two sites 
match up well, and no low-frequency filtering was necessary 
when merging the data. The dominant mode at 10.854 d~' is 
distorted as we detect 2vi and 3vi as well. V2 appears to be 
an independent term, while V3 is equal to v\ + V2- Such asym- 
metries and mode-couplings are common in high-amplitude 
6 Scuti stars. 
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Fig. 5. A single night of data on the BS stars. V6 and V8 were 
observed from La Silla only 



Some of the BS stars fall within the y Doradus instability 
strip as well. V3 is located in the CMD just to the blue side of 
the y Doradus instability strip, but it is nevertheless the only 
one of the six oscillating BS stars which seems to display both 
low- and high-frequency variability. The fast variations in V3 
(Fig. [2] [5] and [6]) are clearly visible. We detect 8 frequencies 
(11-25 d -1 ), with amplitudes in the range 3-56 mmag in B. 
When merging the data from the two sites, it was necessary first 
to remove a relatively strong (25 mmag) signal at 1 d _I from 
the DFOSC B-data only, as mentioned in the caption of Fig. [2] 
The other three data sets (MEROPE BI, DFOSC /) remained 
uncorrected. During the frequency analysis of the merged data, 
we found that three low-frequency terms were present in both 
the B- and /-data. The I/B amplitude ratios of these frequen- 
cies are about 0.4, which is very similar to amplitude ratios of 
the y Doradus stars to be discussed below. This points toward 
oscillations as a possible explanation of the low-frequency vari- 
ability in V3. This result should be confirmed, but it is an in- 
teresting possibility that V3 may oscillate simultaneously in g- 
modes as the y Doradus stars, and in p-modes as the 5 Scuti 
stars. 
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Fig. 7. Amplitude spectra of the new oscillating BS stars V6- 
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4.2. The new variables V6-V8 

In addition to the three oscillating BS stars discussed above, 
we detected short-period oscillations in the three stars V6-V8 
as well. The amplitude spectra are shown in Fig. [7] 
As mentioned above, we only have data from DFOSC on V6. 
As in the case of V3, the B data included a relatively strong 
signal at 1 d -1 , which has been subtracted in Fig. [7] We de- 
tect three frequencies in V6, although v\ does not stand out 
clearly in the / data. There appears to be further frequencies 
present in the range 10-15 d _1 , but these are below the Act de- 
tection threshold of our data. Although V6 is located inside the 
y Doradus instability strip as well, we do not find evidence for 
low-frequency, y Doradus type variability. 

For V7, we detect a single frequency, but the residuals after 
subtracting v\ indicate that further frequencies may be present 
- especially a peak near 9.05 or 10.05 d~' with a B-amplitude 
of about 2 mmag just escapes a Act detection in both the B and 
/-data. V8, on the other hand, is multiperiodic with at least 7 
frequencies, although vj does not meet the Act criterion for de- 
tection. It is, however, present in both the B- and /-filter data 
and is therefore included in Table Q] Further frequencies are 
probably present in the light curves of V8, but escapes detec- 
tion in our data. No low-frequency filtering has been applied to 
the V8-data shown in Fig. [7] 

We finally mention the star V10, which displays variability 
at the mmag-level with a periodicity near 10 d -1 . This star is, 



however, positioned in the CMD near the giant stars of the clus- 
ter and is close to the saturation limit of our data. Little can be 
said at this stage about the nature of this object, the variability 
seems 5 Scuti like, so it is probably not a cluster member. 

4.3. Mode-identification of the oscillating BS stars 

Mode-identification in 5 Scuti stars is a notorious problem. 
Unless the modes can be identified directly, as in the case of 
the double-mode SX Phe stars, or from combined spectroscopy 
and photometry for bright stars such as FG Vir (Viskum et al. 
1998, Breger et al. 1999, Zima et al. 2006), one must turn to 
direct comparison with theoretical models. The latter method 
is further complicated in the present case, as the 6 Scuti vari- 
ables are blue stragglers, for which the formation history is 
unknown (see, e.g., Gilliland et al. 1998; Tian et al. 2006). 
In the case of NGC2506, cluster parameters and isochrones 
will be constrained from the binary analysis of spectroscopic 
(VLT data are already obtained, and will be published in a 
separate paper) and photometric (this paper) data on detached, 
eclipsing cluster binaries. In parallel with this, we will initi- 
ate dedicated model calculations for the oscillating blue strag- 
glers, but this is beyond the scope of the present paper. Instead, 
we will use a period-luminosity relation (PLR) from Petersen 
& Christensen-Dalsgaard (1999, their Eq. 4), along with ex- 
pected period ratios for consecutive, radial overtones in 6 Scuti 
stars, to discuss the detected frequencies (see e.g., Petersen & 
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Table 1. Variability data for the oscillating BS stars. 
Frequencies are in cT 1 and amplitudes are given in mmag. The 
quoted errors on the frequencies are standard deviations on the 
values determined from extensive simulations (see text). 



ID 


Vj 


V 


A, 


A s 


S/N B 


VI 


V\ 


13.630(1) 


3.2 


6.8 


15.1 




V2 


14.684(1) 


2.8 


5.1 


11.2 




V 3 


14.716(1) 


1.7 


2.9 


6.4 




V 4 


16.688(5) 


1.5 


2.8 


6.2 




V5 


16.896(2) 


1.0 


2.2 


5.0 




V6 


17.043(2) 


0.9 


1.8 


4.0 




V 7 


17.672(1) 


2.2 


4.5 


10.0 


V2 


Vl 


10.854(1) 


55.8 


121.4 


215.0 




V2 


20.141(1) 


7.5 


14.2 


29.5 




2v, 


21.709(1) 


8.5 


17.0 


35.3 




v 3 


30.995(2) 


1.2 


2.8 


5.9 




3v, 


32.563(3) 


1.3 


3.1 


6.4 


V3 


Vl 


0.526(1) 


3.2 


6.9 


3.1 




V2 


0.827(6) 


4.4 


10.4 


4.7 




V 3 


0.908(3) 


4.9 


12.2 


5.5 




V 4 


11.243(2) 


1.7 


5.8 


6.4 




V5 


12.031(1) 


7.4 


16.8 


18.6 




Vf, 


12.265(1) 


25.9 


56.2 


62.0 




V 7 


12.748(5) 


2.2 


3.9 


4.4 




V8 


13.088(1) 


10.0 


18.2 


20.1 




v 9 


18.555(4) 


2.5 


3.7 


4.1 




VlO 


24.336(5) 


1.7 


3.2 


4.5 




V\ 1 


24.530(2) 


2.1 


4.3 


6.0 


V6 


Vl 


9.441(3) 


0.9 


3.7 


5.9 




v 2 


10.767(1) 


5.1 


10.9 


17.2 




v 3 


13.958(3) 


1.4 


3.1 


4.9 


V7 


Vl 


10.576(1) 


6.9 


15.3 


22.3 


V8 


Vl 


10.227(2) 


2.9 


4.5 


4.5 




V2 


10.611(1) 


4.4 


11.1 


11.2 




v 3 


10.781(1) 


5.0 


11.7 


11.7 




v 4 


11.550(2) 


4.7 


8.3 


8.4 




V5 


16.149(2) 


3.4 


6.8 


6.8 




V6 


19.423(2) 


3.4 


7.2 


7.2 




v 7 


20.944(3) 


1.8 


3.4 


3.4 



Christensen-Dalsgaard 1996, Gilliland et al. 1998). Although 
the period ratios should be evaluated using dedicated models, 
the ratios are quite constant between models of different mass 
(e.g., Fig. 21 in Gilliland et al. 1998). One should, however, 
be cautious when using period ratios for mode-identifications 
of low-amplitude, multiperiodic 5 Scuti stars, as non-radial 
modes may easily show the same frequency ratios as the ra- 
dial overtones, as noted by Poretti (2003). In addition, one 
should also be cautious of the fact that for the present pa- 
per, we are using period ratios from standard models of non- 
rotating 5 Scuti stars for mode-identification in BS stars. These 
assumptions should be kept in mind when we in the follow- 
ing discuss the mode-identification of the indvidual BS vari- 
ables. In Fig [8] we compare the observed frequencies with the 
period-luminosity relation for the fundamental (F) modes from 
Petersen & Christensen-Dalsgaard (1999) and with consecutive 
radial overtones. 



Kim et al. (2001) identified the dominant frequency of VI 
as the third radial overtone (30), while we find that this fre- 
quency (vi) is better matched with the fourth radial overtone 
(40). We note, by the way, that the frequency ratio between 
vi and V7 is 0.771, which is the expected F/IO frequency ra- 
tio, but such an identification clearly disagrees with the PLR. 
Instead v 7 could possibly be 60, but given the dense frequency 
spectrum, most, if not all, frequencies are likely non-radial. The 
dominant frequency in the high-amplitude 6 Scuti star V2 (vi), 
on the other hand, is in good agreement with a F-mode classi- 
fication, as also noted by Kim et al. The frequency V2 could be 
identified as 30 if we allowed for a 1 d -1 shift to 21.141 d -1 . 
Although we do have significant sidelobes in the amplitude 
spectra, neither the B or / data support such a shift. 

The strongest signal in V3, V6 (at 56 mmag amplitude in E), 
also agrees well with the fundamental mode based on the PLR, 
while all the other frequencies probably are non-radial. Several 
frequencies cluster near 12 d , and three of them (v4,V5,v 7 ) 
form a near equidistant triplet. However, the error simulations 
puts the largest uncertainty on v 7 , of 0.005 d _1 , much smaller 
than the observed deviation from equidistance of 0.071 d" 1 . If 
V3 is a moderate to fast rotator, significant deviations from sim- 
ple triplets are expected (Kjeldsen et al. 1998), but such conclu- 
sions would require comparison with dedicated stellar models 
as well as a spectroscopically derived rotational velocity. 

V6 represents the best case in terms of mode-identification. 
The frequency of the dominant mode (V2) agrees with a F- 
mode classification, and the ratio between V2 and V3 is 0.771, 
which suggests a F/IO identification, vi must then be non- 
radial. V7 is monoperiodic, but the single frequency is in agree- 
ment with the fundamental mode predicted by the PLR, within 
the uncertainties. We do not have Stromgren indices for V8, 
but My ~ 1.15 based on the V-magnitude. The dense spectrum 
near 10 d~' clearly suggests non-radial oscillations, and indeed, 
we have no evidence for radial modes. 

To summarize, comparison with the PLR and expected fre- 
quency ratios suggest radial modes to be excited in VI (40), 
V2 (F), V3 (F), V6 (F/IO) and V7 (F). 

5. Discovery of y Doradus stars 

Among the new variable stars is a population of 15 stars dis- 
playing long-period variability on time scales from hours to 
days. The amplitudes of variability are in the range 5-30 mmag 
and several stars are multiperiodic. They are positioned inside 
the y Doradus instability strip in the CMD (Fig. [3) which, 
along with the periods and amplitudes suggest that they may be 
y Doradus stars (Handler 1999; Handler & Shobbrook 2002; 
Henry et al. 2005). We will in the following argue that these 
stars, labeled VI 1-V25, can in fact be classified as y Doradus 
stars based on the evidence presented below. 

5.1. Evidence for multiperiodic variability 

In Fig. [9] we show light curves for 7 of the 15 y Doradus can- 
didates. They clearly display variability on time scales of days, 
and beating effects resulting from multiperiodicity is evidently 
present in several stars. The amplitude spectra of two of the 
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Fig. 8. Comparison between the observed frequencies in the 
BS stars and a period-luminosity relation for the fundamen- 
tal mode of 5 Scuti stars (thick solid line) and for consecutive 
radial overtones (solid lines) based on theoretical frequency ra- 
tios from stellar oscillation models. The dashed lines indicate 
the lcr error on the PLR and the vertical, dotted lines plots the 
My of the individual variables. The star-symbols indicate the 
frequency of highest amplitude in each star. 



stars are shown in Fig.[10](V17 and V21), along with the spec- 
tral window function shown at the same frequency scale. By 
comparing to the spectral window, which shows the signature 
of a single frequency in the amplitude spectrum, it is again ev- 
ident that several frequencies are present, especially in V17. 
Indeed, we detect 6 and 3 frequencies for V17 and V21, re- 
spectively, in the range 0.1-2.5 d -1 . 

In Fig.QT|we show an example of the merging of data from 
La Palma and La Silla for the variable V19. For 12 of the 15 
stars we have data from both sites and in all cases the overlap- 
ping data agree well (if this was not the case, the star would be 
rejected from the list of variables). This excludes an instrumen- 
tal origin of the long-period variability, which must be intrinsic 
to the stars. For V22,V24 and V25 we have data from La Silla 
only, but the variability is in all cases clearly seen, and present 
in both the B- and /-filter data, leaving little room for doubt 
about the intrinsic nature of the light variability. 

All detected frequencies are listed in Table [2] We find 3-6 
frequencies in several of the stars, and only three stars appear 
mono-periodic in our data. In general we find I/B amplitude ra- 
tios of 0.4-0.6 (Aj/Ab) while a wide range in phase differences 
are observed. The typical residual noise level in the amplitude 
spectra in the frequency range - 6 d -1 are 0.5 - 2.0 mmag for 
these V = 15-16 stars. These numbers explain why Kim et al. 
(2001) could not find y Doradus stars in their dataset. Indeed, 
their dataset consisted of 304 images, and they quote a detec- 
tion limit of 10 mmag at V = 14. The corresponding number for 
our dataset is 2 mmag, as residual noise levels in the amplitude 
spectra of the V=14 stars in our sample are typically below 0.5 
mmag. 
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Fig. 9. Light curves for 7 of the long-period variables 
(y Doradus) stars in our sample. 

5.2. Are theyy Doradus stars? 

To be considered a bona fide y Doradus star, a star should be 
located in the CMD inside the y Doradus instability strip, be 
variable on time scales from ~ 0.4 - 3 days, and all mecha- 
nisms for variability other than oscillations must be excluded. 
Such other mechanisms would typically be binarity (ellipsoidal 
variability) or rotational modulation (star spots). 

All 15 stars are positioned inside the y Doradus instability 
strip near the cluster main sequence turn-off. From isochrone 
fitting we find that the bluemost point on the cluster sequence 
corresponds to a mass of 1.4 - 1.5M . This result is corrobo- 
rated by a preliminary analysis of the binary V4, which places 
the primary at the top of the cluster main sequence with a mass 
of 1.5 + 0.1 M and the secondary further down the main se- 
quence, below the instability strip, with a mass of 1.2 + 0.1 M 
(see Fig. [TBI . This indicates that the turn-off stars are of early 
f -typqj, in good agreement with the y Doradus phenomenon 
occurring in stars of type Al - F5 (Kaye et al. 1999). The po- 
sition of the 15 variables near the cluster turn-off is thus fully 
in agreement with a y Doradus classification. 

1 Allen's Astrophysical Quantities, 4th edition, Arthur N. Cox, Ed. 
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Fig. 11. Three nights of merged data for the y Doradus star 
V19. Filled circles are data from La Palma, open circles from 
La Silla. The amplitude difference between B and / is evident. 



For the multiperiodic variables in our sample, mechanisms 
other than oscillations can readily be excluded as the cause of 
variability (see Handler 1999, Handler & Shobbrook 2002 for 
a detailed discussion). In short, ellipsoidal variables and rota- 
tionally modulated stars would be characterized by a single fre- 
quency, or by two frequencies with a harmonic relation. In the 
case of spots, two close frequencies may appear. For ellipsoidal 



variables, amplitude ratios between different passbands are ex- 
pected to be unity, while V/B amplitude ratios of 0.88-0.89 are 
expected for starspot variables (Henry et al. 2005). For oscilla- 
tions, this ratio is expected to be 0.67-0.80 (Garrido 2000) and 
typically 0.3-0.6 for I/B (see Sect. 5.2), which is also what we 
find for most frequencies (Table|2]l. 

Among the 15 y Doradus stars, there are three mono- 
periodic stars, but for both V13 and V24, further frequencies 
are likely present and only just escape detection in our data set. 
For V13 and V24, the amplitude ratios are also near 0.4, point- 
ing towards oscillations and clearly rejecting ellipsoidal vari- 
ability. The near detection of further frequencies, along with 
the low I/B amplitude ratios suggest that also V13 and V24 are 
oscillating. For V12, the amplitude ratio is 0.7 which rejects el- 
lipsoidal variability. V13 (along with VI 2) are in the Mercator 
field and the data from March and April show that the variabil- 
ity is coherent over time scales of many months. Star spots are 
not expected in y Doradus stars due to the early F-type clas- 
sification (Kaye et al. 1999, Henry et al. 2005), except if they 
are chemically peculiar. As we will show below, none of the 
three monoperiodic stars appear chemically different from the 
rest of the group based on the Stromgren photometry, but all 
15 stars have low values of Sc\, which is an indicator of chem- 
ical peculiarity (Joshi et al. 2006). Thus, we cannot completely 
rule out that V12 could be a rotationally modulated Ap star, 
this can only be done spectroscopically. For V22, for which 
we do not have Stromgren indices except for the position in 
the CMD inside the instability strip, vi and V2 are quite close 
to be harmonically related, but the I/B amplitude ratios of 0.4 
excludes ellipsoidal variability. The evidence therefore suggest 
oscillations as the cause of variability for 14 of the 15 stars, in 
agreement with a y Doradus classification. 

5.3. Comparison with the known y Doradus stars 

We have compared the Stromgren indices of the 15 NGC2506 
variables to those of the bona fide and prime-candidate 
y Doradus stars based on updated tables kindly provided by G. 
Handler. For the latter stars we extracted Stromgren data from 
the catalogue of E. H. Olsen (private communucation) and used 
the b-y calibration of Crawford (1975) to derive Sc\ and Snt\. 
The result is shown in Fig. [T2l 

In the My,{b - y)o diagram, the NGC2506 stars overlap 
completely with the known variables. In the (Sc\,6m{) diagram, 
the NGC2506 stars have a slightly more negative 6c\ and, in 
general, slightly higher 8m\. The higher 5m \ values suggest 
lower metallicity, in agreement with the metal content of the 
cluster, while the more negative 6c \ suggest that the stars may 
be more chemically peculiar and therefore could be Ap stars 
(Joshi et al. 2006). The NGC2506 stars are, however, still lo- 
cated very close to the known y Doradus stars in this (Sc\,Sm\) 
diagram, where we have marked the variables with three or 
more frequencies with an open circle as well. The variability 
of those multiperiodic stars are certainly due to oscillations, 
and this subset includes the two stars with the most negative 
values of Sc\. The 15 NGC2506 stars all have very similar in- 
dices; they form a homogeneous group in both the CMD and 
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the {5c\,5m\) diagram, and nothing suggests that any one star 
is chemically different from the others. We mention in pass- 
ing that this supports oscillations as the cause of variability in 
V12 as well, as the star has indices similar to those of the 14 
other y Doradus candidates. The comparison with the known 
oscillators therefore strengthen the y Doradus classification of 
the NGC2506 stars. The minor shift of the NGC2506 stars in 
5c\ ,5mi may be due to a small zero-point shift of our Stromgren 
photometry compared to the known variables, which are typi- 
cally closer field-stars, or to an actual increase in the 5c\,5m\ 
parameter space of the y Doradus stars in general. 

The fact that all the NGC2506 y Doradus stars lie on or near 
the cluster main sequence inside the y Doradus instability strip 
is a strong indication that they are all cluster members. In fact, 
plots of various combinations of the Stromgren indices given in 
Table [3] all suggest that the y Doradus stars are cluster member 
turn-off stars, in full agreement with the CMD in Fig. [3] If some 
of the 15 stars were foreground or background objects, and thus 
were of a different spectral type, they would have Stromgren in- 
dices separating them out from the large group of cluster turn- 
off stars. It would also be difficult to explain the pulsational 
behaviour of such objects; background Cepheids or RR Lyrae 
stars would have a characteristic light curve with a larger ampli- 
tude of variability. A reddened, background SPB variable could 
show the same type of variability as in the NGC2506 stars, but 
would be unlikely to show Stromgren indices identical to those 
of F-stars in a metal-poor cluster. 

Furthermore, we have VLT/FLAMES spectra for V20, V21 
and V24, of which the first two are multiperiodic oscillators, 
in a data set of 120 stars in the cluster. We have derived pre- 
liminary radial velocities based on these spectra (see Sect. 6), 
and the result is that at least V20 and V24 are cluster members 
as their mean radial velocity corresponds to the cluster mean 
velocity. This provides additional, strong support for the clas- 
sification of all 15 stars as cluster member 7 Doradus stars. The 
velocities of V21 are in agreement with cluster membership but 
long-term RV variations may be present. 

We therefore suggest that the NGC2506 variables VI 1- 
V25 except V12 should be included in the list of bona fide 
y Doradus stars. V12 should be regarded as a prime candidate 
until the nature of the variability has been established with cer- 
tainty. This population of y Doradus stars makes NGC2506 ex- 
tremely interesting for studying the y Doradus phenomenon, as 
these stars cover a large fraction of the instability strip and at 
the same time are located very near to the evolutionary point of 
turn-off from the main sequence. 

We finally note that we did not find 6 Scuti oscillations in 
any of the y Doradus stars, despite the fact that many of them 
are inside the 6 Scuti instability strip as well. This supports the 
idea that y Doradus and S Scuti oscillations are rarely present in 
the same star (Handler & Shobbrook 2002, Henry et al. 2005). 
Both types of oscillation have, however, been found in HD8801 
(Henry & Fekel 2005) and in BD+184914 (Rowe et al. 2007) - 
and possibly in V3 in this paper. It is interesting to note that we 
do not find S Scuti oscillations in any of the turn-off stars, de- 
spite many candidates inside the 6 Scuti instability strip. In con- 
trast, the open cluster NGC1817 also has turn-off inside (in the 
middle of) the 6 Scuti instability strip, but that cluster, which is 
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Fig. 12. Comparison of stellar parameters between the 
NGC2506 y Doradus stars and the known variables (Handler, 
private communication). The bona fide y Doradus stars are 
shown as squares, the prime candidates as triangles, the 
NGC2506 stars as filled circles, and the instability strip is from 
Handler & Shobbrook (2002). In the bottom panel, the mul- 
tiperiodic (3 or more frequencies) NGC2506 stars have been 
marked with an open circle as well. 



half as old as NGC2506, is hosting a population of at least 12 
turn-off 6 Scuti stars (Arentoft et al. 2005). 

5.4. Comparison with theoretical predictions 

As mentioned in the introduction, y Doradus stars are rich stel- 
lar laboratories to study astrophysical phenomena. This, how- 
ever, makes them complex and challenging targets for which 
we need as many additional observational constraints as possi- 
ble. Such constraints will be obtained when we determine with 
high precision the masses and radii of the two eclipsing bina- 
ries in NGC2506. The high-resolution spectroscopic datasets 
involved, obtained with the VLT at Paranal Observatory, are 
currently being reduced and analysed. The datasets will pin- 
point the cluster isochrone with an unprecedented precision for 
NGC2506, which will lead to masses, radii and the age of all 
stars along the cluster main sequence, including the y Doradus 
stars. In addition, a detailed abundance analysis of cluster stars 
based on VLT/FLAMES spectra will shed light on the metallic - 
ity of the cluster stars. These results, which will be published in 
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Table 2. Variability data for the y Doradus stars. Frequencies 
are in d -1 , amplitudes in mmag and phases in degrees. The 
quoted errors on the frequencies, amplitude ratios and phase 
differences are standard deviations on the values determined 
from extensive simulations (see text). 



ID 


Vj 


V 




S/N B 


A ; /A B 


</>, - (f> B 


Vll 


Vl 


1.165(1) 


7.93 


6.78 


0.41(4) 


-1(5) 




Vl 


1.270(1) 


5.12 


4.37 


0.40(6) 


6(8) 




V\ 


1.400(2) 


8.58 


7.34 


0.40(3) 


3(5) 


V12 


Vl 


1.395(4) 


9.26 


5.84 


0.70(6) 


30(6) 


V13 


V\ 


0.830(1) 


8.10 


6.00 


0.39(4) 


13(6) 


V14 


V] 


0.758(4) 


2.130 


4.57 


0.68(13) 


33(12) 




V2 


0.870(1) 


17.750 


38.08 


0.43(2) 


-8(2) 




V3 


0.927(1) 


4.090 


8.77 


0.31(8) 


88(18) 


V15 


V] 


0.495(4) 


3.60 


3.16 


0.46(13) 


35(17) 




Vl 


0.614(1) 


14.93 


13.10 


0.48(2) 


4(2) 




v 3 


1.130(1) 


6.28 


5.51 


0.45(4) 


131(5) 


V16 


V] 


0.474(20) 


9.599 


5.02 


0.29(3) 


14(13) 




Vl 


0.795(1) 


20.603 


10.78 


0.52(2) 


-13(2) 




V3 


0.918(3) 


11.459 


6.00 


0.49(4) 


1(4) 




v 4 


1.730(13) 


6.144 


3.21 


0.73(11) 


-4(8) 


V17 


V] 


0.136(1) 


10.875 


7.64 


0.68(4) 


21(3) 




Vl 


0.606(1) 


28.233 


19.85 


0.46(2) 


4(2) 




V3 


0.645(7) 


15.658 


11.01 


0.51(3) 


-11(3) 




V4 


1.268(1) 


7.362 


5.17 


0.55(5) 


-3(5) 




V5 


1.537(2) 


7.213 


5.07 


0.35(11) 


28(15) 




V6 


2.439(1) 


5.711 


4.01 


0.65(8) 


-2(7) 


V18 


V\ 


0.206(1) 


8.939 


9.51 


0.80(4) 


-15(3) 




vz 


0.674(2) 


3.767 


4.01 


0.96(11) 


-14(6) 


V19 


V] 


0.156(1) 


7.190 


6.13 


0.56(4) 


1(4) 




Vl 


0.442(2) 


15.774 


13.44 


0.51(2) 


1(2) 




V3 


0.583(1) 


20.767 


17.70 


0.39(2) 


0(3) 




V4 


1.344(3) 


6.928 


5.90 


0.51(5) 


-27(5) 




V5 


1.593(5) 


6.902 


5.88 


0.51(6) 


5(6) 


V20 


V\ 


0.842(25) 


7.610 


11.23 


0.43(5) 


-28(7) 




Vl 


0.883(1) 


9.963 


14.70 


0.66(5) 


-2(3) 




V3 


1.701(16) 


3.011 


4.44 


0.47(8) 


-28(11) 


V21 


Vl 


0.503(1) 


5.069 


3.50 


0.53(9) 


4(15) 




v 3 


0.648(1) 


17.793 


12.30 


0.40(2) 


-15(2) 




v 3 


1.740(1) 


5.434 


3.76 


0.60(5) 


29(6) 


V22 


V\ 


0.700(1) 


14.214 


6.18 


0.44(3) 


-8(4) 




Vl 


1.378(1) 


12.684 


5.52 


0.42(3) 


2(4) 


V23 


V\ 


1.168(3) 


9.616 


6.35 


0.47(4) 


1(5) 




Vl 


1.260(3) 


14.415 


9.52 


0.52(3) 


12(3) 


V24 


V\ 


0.916(1) 


14.458 


12.25 


0.45(3) 


16(5) 


V25 


V\ 


1.178(1) 


14.747 


10.84 


0.60(3) 


12(2) 




Vl 


1.987(2) 


7.055 


5.19 


0.50(5) 


-18(8) 



a separate paper, will provide the needed, additional constraints 
and we therefore postpone a detailed seismic analysis of the y 
Doradus candidates until they are available. 

We did, however, perform a preliminary theoretical model- 
ing to verify whether there is at least a crude correspondence 
between the observed amplitude ratios and phase differences 
and the theoretical ones. Rather than modeling each target in- 
dividually and being subjected to the possible large, systematic 
uncertainties on their positions in the HR-diagram, we choose 
to compute 4 generic stellar models using the stellar evolution 
code CLeS (R. Scuflaire, University of Liege, private commu- 







2.0 M 












" .1.8 Mq 






.1.6 M '- 












1.4 Mq r 




■ / 











3.9 3.88 3.86 3.84 3.82 3.8 3.78 3.76 3.74 



Log(T eff ) 

Fig. 13. An HR-diagram showing the 4 generic y Doradus mod- 
els as squares together with the positions of the y Doradus can- 
didates as crosses. Also shown are the ZAMS, evolutionary 
tracks and the blue and red border of the y Doradus instabil- 
ity strip for { = 1 modes and for the mixing length parameter 
a = 2. 

nication). Each of these models was computed for the same 
cluster age as specified in Section [2] but with different effec- 
tive temperatures going from the blue to the red border of the y 
Doradus instability strip. These generic models, together with 
the positions of the y Doradus candidates are shown in the HR- 
diagram in Fig. [T3l 

Consequently, we computed the i = 1,2 and 3 eigen- 
frequency spectra for each of the generic models using 
the non-adiabatic oscillation code MAD. This code uses 
time-dependent convection (TDC) theory as described by 
Grigahcene et al. (2005) and Dupret et al. (2005a), and was 
successfully applied to five y Doradus field stars by Dupret et 
al. (2005b). From these spectra, and for each degree t, we se- 
lected the modes with eigenfrequencies closest to the observed 
frequencies and computed the non-adiabatic relative amplitude 
and phase of the T e g and gravity variation. These quantities 
were then used to compute the photometric amplitude ratio I/B 
and phase difference I - B. The theoretical model used is the 
same as the one used in Dupret et al. (2005b), and incorpo- 
rates Kurucz (1993) atmosphere models and the 4-coefficient 
non-linear limbdarkening law of Claret (2000). Through inter- 
polation in metallicity grids, the results were computed for the 
cluster metallicity mentioned in Section|2]([Fe/H] = -0.20). 

Three important results come out from our comparison 
between theoretical predictions and observations. First, the 
generic models are unable to explain the observed oscillation 
frequencies below 0.25 d . In particular, v\ = 0.136 d -1 of 
V17, vi = 0.206 d 1 of V18, and vj = 0.156 d~' of V19 
are predicted to be stable (i.e. not excited). The most impor- 
tant parameter determining the (in)stability of eigenmodes in 
y Doradus stars is the mixing-length parameter a. We used 
the value a = 2.0 because a) testing revealed that this value 
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leads to the largest number of unstable modes, and b) Dupret 
et al. (2005b) successfully used this value to explain best the 
observations of five field y Doradus stars. But in the case of 
our cluster targets this value seems not able to explain the com- 
plete observed frequency spectrum. We remark, however, that 
the observed frequency V3 = 0.343 d _1 of the well known y 
Doradus star 9 Aur is also predicted to be stable by Dupret et 
al. (2005b). 

The second result concerns the photometric phase differ- 
ences I - B. Fig. [14] shows the theoretically predicted photo- 
metric amplitude ratios I/B and corresponding phase difference 
I — B for all modes of all targets, assuming degrees £=1,2 and 
3. Comparing this figure with the corresponding observational 
Fig- reveals a striking difference: for many of the oscillation 
modes we observe a positive phase difference I — B which can- 
not be explained by the models. The phase difference of these 
modes is usually positive well within the error bars, and this is 
even so for several of the highest-amplitude modes for which 
no aliasing-problems was detected (red triangles). There is lit- 
tle to compare to in the literature but we mention the relevant 
information in the forthcoming paper of Cuypers et al. (2007, in 
preparation). In this paper a sample of 19 y Doradus field stars 
is presented that was photometrically monitored in the Geneva 
passbands. The authors always find a phase difference between 
the U and G filter that is very close to zero (within the stan- 
dard errors). This turns out to be in agreement with the theory. 
Indeed, we performed similar simulations for the Geneva U 
and G filters, and we found that the phase difference U - G be- 
haves somewhat different than I - B. First, it can theoretically 
be both positive and negative. Secondly, regardless the degree, 
the phase difference U - G is almost always close to zero. The 
results of Cuypers et al. are therefore not surprising. It does 
seem to show, however, that the phase difference / - B is more 
suitable for mode identification as it has more discriminating 
power. 

A third remark worth mentioning is the scatter of the ob- 
served modes in Fig. [4] Intuitively we would expect more £ = 1 
modes to be present among our sample as a selection effect be- 
cause these modes suffer less surface cancellation. This effect 
seems to be present in the sample of Dupret et al. (2005b) who 
identified 8 out of 12 studied y Doradus modes unambiguously 
as dipole modes. Note, however, that the I — 1 modes in Fig. [14] 
(blue circles) occupy a remarkably small region given the fact 
that the 4 generic models from which they were computed are 
spread over a significant part of the y Doradus instability strip. 
The rather large scatter in the observational diagram in Fig. [4] 
therefore shows that the majority of our sample of y Doradus 
modes cannot easily be identified as dipole modes. 

Given the evidence, we should conclude that either the 
observational uncertainties on the photometric amplitude ra- 
tios and phase differences are significantly larger than we es- 
timated, or the theoretical model is inadequate. We verified 
that changing the metallicity (e.g. to [Fe/H] = 0.0 or [Fe/H] 
= -0.3) does not change qualitatively our results. In our forth- 
coming theoretical study we will explore a larger parameter 
space where we will examine in more detail the influence of 
the overshoot parameter, and the rotational velocity. The latter 
in particular may be important for y Doradus stars as for these 
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Fig. 14. The theoretical photometric amplitude ratios I/B ver- 
sus phase difference / - B. It concerns the modes of the 4 
generic models shown in Fig.[l3]for which the eigenfrequency 
is close to one of the observed frequencies in Table|2] The blue 
circles denote t — 1 modes, the red squares I — 2 modes and 
the yellow triangles I — 3 modes. The results were computed 
for the cluster metallicity stated in Section|2] and for the mixing 
length parameter a — 2. 

stars the rotation over oscillation frequency ratio can easily be 
comparable to or larger than 1 . 

6. Eclipsing Binaries 

As mentioned above, NGC2506 hosts the two detached, eclips- 
ing binaries V4 and V5. We have used the present photometric 
data set to determine, for the first time, the periods of these two 
binaries: 2.868 and 10.078 days for V4 and V5, respectively. 
The phased light curves are shown in Fig. [15] including some 
nights of poor quality. The B light-curve of V4 display minor 
offsets between data taken at the same phase, but on different 
nights. The same is the case for /, although to a lesser extent. 
Fig.[T6lillustrates this in more detail. If we assume that the bi- 
nary is a cluster member, as indicated by the spectroscopic data 
to be discussed below, the primary component should lie on the 
cluster sequence near turn-off, i.e., somewhere along the upper- 
most of the green curves in Fig. [16] lower panel. If we also re- 
quire the secondary to be on the cluster sequence, and that the 
combined light of the two main sequence stars should be equal 
to the observed, total light, we obtain a photometric solution 
that puts the primary inside the y Doradus instability strip and 
the secondary somewhat below it, as indicated by the red dots 
in the diagram. In the upper panel of Fig. [16] we show the light 
curve in more detail. Additional variability is present on a time 
scale of several days which is, however, longer than the typical 
y Doradus type variability. Further discussion of this object is 
therefore postponed until all data, including the spectroscopy, 
are available. 

From the photometric data, we also detect four new eclips- 
ing binaries (V9 and V26-V28). Light curves are shown in 
Figs. [17] [18] and [19] and we refer to Table [3] for photomet- 
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nc indices. V26 is a giant star, near the limit of saturation in 
our data, especially in I. Fig. [17] plots the evidence of bina- 
rity. V27 (Fig. |T8T> is possibly a new, detached eclipsing binary 
positioned in the CMD on the cluster main sequence. A single, 
probable eclipse was observed in both B and / during one night. 
Finally, we detect the two binaries V9 and V28 (Fig. [19]). V9 is 
located in the CMD inside the y Doradus instability strip, near 
the cluster sequence, and may be an Algol- or /3 Lyrae type 
binary. V28 is faint, red and, from the position in the CMD 
far from the cluster main sequence, not a cluster member. The 
short period (0.2199 days), and the light-curve shape, suggest 
a W UMa classification. 

The spectroscopic data-set mentioned above consists of 15 
epochs of VLT UVES spectra for 5 stars (V4 and four giants) 
and Giraffe/FLAMES spectra of 120 stars (including V5) ob- 
served in service mode. The aims of these data are to anal- 
yse the binary stars for deriving precise, absolute parameters, 
as well as to use the Giraffe spectra to determine metallic- 
ity/abundances, ages and cluster membership. The latter allows 
a cleaning of the bright part of the CMD for non-members, 
which in turn allows much stronger constraints to be put on the 
isochrones. The analysis of the NGC2506 spectra is still ongo- 
ing and will be the subject of a later paper, but we mention that 
a preliminary analysis of the UVES spectra shows that V4 is a 
double-lined binary and a cluster member. 

We have also done a preliminary reduction of the Giraffe 
spectra to derive approximate radial velocities for the 120 ob- 
served stars, among these V5 and V9. Both stars show clear 
velocity variations centered on the cluster mean value of ap- 
proximately +84 km/s. The two y Doradus stars V20 and V24, 
for which we have Giraffe spectra, have constant velocities over 
the 15 epochs, with values near the cluster mean. They are thus 
cluster members, as are V4, V5 and V9. For V2 1 , long-term RV 
variations are possibly present, with values between +70 and 
+90 km/s, but further analysis is needed before any firm con- 
clusions can be drawn about the membership status and possi- 
ble dublicity of this object. 

7. Summary and conclusions 
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Fig. 15. Phased light curves for the detached binaries V4 and 
V5. 
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Our two-colour, dual-site photometric campaign from La 
Palma and La Silla was a success as we have quadrupled the 
number of known variables in the open cluster NGC2506. 
Based on our results we conclude that the cluster contains (at 
least) 6 oscillating BS stars, 15 y Doradus stars, and 6 eclipsing 
binaries. 

Prior to this work, only three short-period variables were 
known in NGC2506, and only one of them was classified as 
a oscillating BS star. We now know 6 oscillating BS stars, of 
which 5 are multi-periodic. One BS star (V3) shows both short - 
and long-period variability. To our knowledge, this would be 
the first BS star for which both p-mode and g-mode oscilla- 
tions are discovered. Given the fact that we detected 1 1 oscil- 
lations frequencies in this star, makes this blue straggler a very 
promising target for follow-up research. We therefore reached 
our goal of data gathering and selection of interesting targets 
suitable for future investigation. We finally mention that a pre- 
liminary mode identification seems to indicate that in four of 
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b-y 

Fig. 16. 24-day light curve (top) and photometric solutions to 
the V4 binary system. The solid red lines are the y Dor insta- 
bility strip. The cross denotes the position of V4 in the CMD 
and the green curves are possible photometric solutions of the 
two components. The red dots mark a solution that puts both 
components on the cluster sequence. 

the BS stars the fundamental radial mode is excited, and for one 
star ( V6) there is evidence of the fundamental and first overtone 
frequencies being excited. 

We report the discovery of 15 y Doradus stars of which 
most are multiperiodic, and at least 14 of them are bona fide 
y Doradus stars, while one is a prime y Doradus candidate. The 
(b-y) colours, My and the 6c\,5m\ indices agree well with the 
corresponding values for the known y Doradus variables. 
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Table 3. Stellar parameters for the variables in NGC2506. Type names refer to oscillating blue straggler for BS, detached 
eclipsing binary for DEB, contact binary 

for CB, y Doradus for GD, and eclipsing binary (unspecified) for EB. V10 has not been classified, and Stromgren photometry is 
not available for V8, V22 and V28 . 



ID 


<?2000 


^2000 


V 


b-y 


v - y 


m, 


Cl 


(b-yh 


Co 


T* 


My 


Mm 


BC 


logL 


Type 


VI 


7 59 58.1 


-10 45 54 


13.660 


0.150 


0.468 


0.168 


0.994 


0.096 


0.983 


8330 


0.965 


0.903 


-0.062 


1.537 


BS 


V2 


7 59 53.6 


-10 45 48 


14.510 


0.185 


0.548 


0.178 


0.819 


0.131 


0.808 


7950 


1.815 


1.783 


-0.031 


1.185 


BS 


V3 


7 59 57.9 


-10 47 20 


14.727 


0.207 


0.561 


0.147 


0.820 


0.153 


0.809 


7740 


2.032 


2.009 


-0.022 


1.095 


BS 


V4 


8 00 08.2 


-10 45 50 


14.645 


0.292 


0.706 


0.122 


0.586 


0.238 


0.575 


6920 


1.950 


1.900 


-0.050 


1.138 


DEB 


V5 


8 00 10.3 


-10 43 17 


17.430 


0.456 


1.133 


0.221 


0.199 


0.402 


0.188 












DEB 


V6 


8 00 17.2 


-10 45 35 


14.677 


0.243 


0.701 


0.215 


0.651 


0.189 


0.640 


7360 


1.982 


1.965 


-0.017 


1.113 


BS 


V7 


7 59 53.4 


-10 47 12 


14.694 


0.220 


0.599 


0.159 


0.752 


0.166 


0.741 


7600 


1.999 


1.983 


-0.015 


1.105 


BS 


V8 


8 00 34.5 


-10 46 50 


13.80 


0.22 






















BS 


V9 


8 0005.4 


-10 43 39 


15.289 


0.268 


0.666 


0.130 


0.681 


0.214 


0.670 


7160 


2.594 


2.563 


-0.030 


0.873 


EB 


V10 


7 59 54.9 


-10 46 46 


13.079 


0.522 


1.287 


0.243 


0.513 


0.468 


0.502 












? 


VII 


8 00 04.6 


-10 48 43 


15.454 


0.275 


0.681 


0.131 


0.621 


0.221 


0.610 


7080 


2.759 


2.720 


-0.038 


0.810 


GD 


V12 


8 00 05.0 


-10 47 11 


15.508 


0.269 


0.654 


0.116 


0.654 


0.215 


0.643 


7140 


2.813 


2.780 


-0.032 


0.786 


GD 


V13 


8 00 04.7 


-10 46 39 


15.599 


0.269 


0.663 


0.125 


0.621 


0.215 


0.610 


7130 


2.904 


2.870 


-0.034 


0.750 


GD 


V14 


7 59 55.8 


-10 48 41 


15.376 


0.262 


0.687 


0.163 


0.610 


0.208 


0.599 


7180 


2.681 


2.652 


-0.029 


0.838 


GD 


V15 


7 59 56.5 


-10 47 03 


14.983 


0.267 


0.676 


0.142 


0.566 


0.213 


0.555 


7120 


2.288 


2.253 


-0.034 


0.997 


GD 


V16 


7 59 57.7 


-10 46 09 


15.333 


0.243 


0.624 


0.138 


0.668 


0.189 


0.657 


7370 


2.638 


2.621 


-0.017 


0.850 


GD 


V17 


8 00 07.4 


-10 43 33 


15.521 


0.257 


0.662 


0.148 


0.621 


0.203 


0.610 


7230 


2.826 


2.802 


-0.024 


0.778 


GD 


V18 


8 00 01.0 


-10 45 06 


15.517 


0.254 


0.642 
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Fig. 17. The new probable binary V26. The red curve shows a 
2-frequency fit with a 1:2 relationship between them. 



We have compared frequencies, amplitude ratios and phase 
differences of the y Doradus stars to state-of-the-art y Doradus 
models. The results of this comparison are that 1) the lowest, 
detected frequencies (below 0.25 d _1 ) are predicted to be stable 
in the models, 2) we observe positive I — B phase differences 
in disagreement with the models, and 3) the majority of the 
observed frequencies appear not to be dipole modes, contrary 
with what is expected, as cancellation effects are stronger for 
modes of higher degree. These results clearly ask for further 
examination. We plan a forthcoming theoretical study and ad- 
ditional spectroscopic data to investigate the effect of stellar 
rotation on the observables. 

We have for the first time determined periods for the two 
detached, eclipsing binaries V4 and V5, and preliminary spec- 
troscopic results confirms cluster membership for both. 





-0 


05 


o 





00 







05 







10 




-0 


05 


o 





00 













05 







1 



1 V27 




** 


% i : 

° 


— i i i i i 


i''' 


i ' ' 

• i 


1 



























HJD-53375 



Fig. 18. The new probable binary V27; four nights of data are 
shown, with a possible eclipse taking place during the last 
night. 



All these results combined strongly confirm and increase 
our interest in NGC2506 as an excellent target cluster for study- 
ing BS and y Doradus stars. 
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The latter object is not a cluster member based on the position 
in the CMD. 
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Fig. 20. DFOSC image serving as a finding chart for NGC2506, with the variables indicated. The image orientation is indicated 
by arrows. 



